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Abstract 

This paper discusses the ATLAS potential to study Supersymmetry in the "focus point" 
region of the parameter space of mSUGRA models. The potential to discover a deviation 
from Standard Model expectations with the first few fb 1 of LHC data was studied using the 
, , parametrized simulation of the ATLAS detector. Several signatures were considered, involving 

hard jets, large missing energy, and either &-tagged jets, opposite-sign isolated electron or muon 
pairs, or top quarks reconstructed exploiting their fully hadronic decays. With only 1 fb -1 of 
data each of these signatures may allow observing an excess of events over Standard Model 
expectation with a statistical significance exceeding 5 standard deviations. Furthermore, each 
of the two invariant mass distributions of the two leptons produced by the \3 ~* Xi^ + '~ an d 
the X2 — * Xi^ + '~ three-body decays has a kinematic endpoint which measures the difference 
between the masses of the parent and daughter neutralino. An analytical expression was 
, derived for the shape of this distribution and was used to fit the simulated LHC data. A 

O • measurement of the x® — X? an d X3 — Xi mass differences was obtained and this information 

' was used to constrain the MSSM parameter space. 

1 Introduction 

> 

l/*-) ' One of the best motivated extensions of the Standard Model is the Minimal SuperSymmetric 

Model (MSSM) pQ. Because of the large number of free parameters of the general MSSM, the 
studies in preparation for the analysis of LHC data are mostly performed in a more constrained 
framework, obtained making some assumptions about the breaking mechanism of Supersymmetry. 
Most studies are performed in the minimal SUGRA framework [I], which has five free parameters: 
the common mass mo of scalar particles at the grand-unification energy scale, the common fermion 
mass m 1 / 2 , the common trilinear coupling Ao, the sign of the Higgsino mass parameter \i and the 
ratio tan/3 between the vacuum expectation values of the two Higgs doublets. 

Constraints are provided by searches made by experiments at accelerators (in particular LEP [5]) 
and by the requirements that the radiative electroweak symmetry breaking is consistent with the 
Standard Model. A strong point of Supersymmetry is that in case of exact R-parity conservation 
the lightest SUSY particle (LSP) is stable and can thus provide a candidate for Dark Matter. 
Because of cosmological considerations the LSP must be neutral and weakly interacting and in 
mSUGRA the suitable candidate is the lightest neutralino Xv It is therefore natural to apply 
the additional constraint that the neutralino relic density fi^ in the present universe should be 
compatible with the density of non- baryonic Dark Matter, which is n DM h 2 = 0.105^°;°^ [3J. If 
there are other contributions to the Dark Matter one may have £1% < Qdm- 

In most of the mSUGRA parameter space, however, the neutralino relic density is larger than 
dm 0] ■ An acceptable value of relic density is obtained only in particular regions of the parameter 
space, noticeably: 

• in a region with a relatively low value of the SUSY mass scale [bulk region). 

• for m 1 / 2 >> wio, when the mass of the scalar r is close to the mass of the lightest neutralino, 
so that xt annihilation in the early universe reduces the relic density (f co- annihilation 
region) . 

• for large value of tan /3, there is a funnel in the parameter space where the mass of the pseudo- 
scalar Higgs boson is nearly twice the one of the neutralino, enhancing the XX annihilation 
cross section (Higgs funnel region). 
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• for mi/2 << mo a region exists [5], where the lightest neutralino has a significant Higgsino 
component, enhancing the XX annihilation cross section (focus point region). 

In this paper a study of the ATLAS potential to discover and study Supersymmetry in the 
focus point scenario is presented. 

The paper is organized as follows. In Section [2] a scan of the minimal SUGRA parameter 
space is performed to map the focus point region with an acceptable relic density. The theoretical 
uncertainties on the SUSY mass spectrum and the variation of the masses across the focus point 
parameter space are discussed. Finally, a benchmark point is selected for more detailed studies. 
For this benchmark, the production of Supersymmetric particles is dominated by the electroweak 
production of neutralinos and charginos, which is not easily observable in a hadronic collider, 
and by the pair production of gluinos, followed by the decay of each gluino into (mostly) third- 
generation quarks and a neutralino or chargino. The gluino pair production results in events 
with hard jets, including 4 6-jets per event, large missing energy from the lightest neutralinos 
which escape detection, and possibly isolated leptons from neutralino and chargino decays. These 
signatures are used in the studies reported in Section [3] to 

In Section [3] the performance of inclusive search strategies based on the presence of hard jets, 
large missing energy, and 6-tagged jets to discriminate the SUSY signal from the Standard Model 
background is studied. 

In Section [4] the reconstruction of the kinematic edge of the invariant mass distribution of the 
two leptons from the decay Xn ~ * Xi' + ^~ is discussed. The presence of isolated lepton pairs (in 
addition to hard jets and missing energy) is a promising discovery channel with the first few fb _1 
of LHC data; with a larger dataset, the kinematic endpoints of the two-lepton invariant mass dis- 
tribution allow the measurement of two constraints on the masses of the three lightest neutralinos. 
To this purpose, an analytical expression for the shape of the invariant mass distribution of the 
lepton pairs arising from the three-body leptonic decay of the neutralinos has been derived, in the 
hypothesis of heavy scalar masses (which is a suitable approximation for the focus point). The 
resulting function was used in the fit of the dilepton invariant mass distribution obtained with 
simulated data, obtaining a measurement of the X2 ~ Xi an d X3 — Xi mass differences. 

In Section[S]the reconstruction of the gluino decays is discussed. To this purpose, the top quark 
fully-hadronic decays are reconstructed; the presence of an excess of tb and tt pairs is a possible 
discovery channel while their invariant mass places constraints on the gluino mass scale. 

Finally, Section [(^investigates which constraints one can put on the Supersymmetry parameters 
from the measurements of the two neutralino mass differences. 

2 Scans of mSUGRA parameter space 

In order to find the regions of the mSUGRA parameter space which have a relic density compatible 
with cosmological measurements, the neutralino relic density was computed with micrOMEGAs 
1.31 [BJ, interfaced with either ISAJET 7.71 or SOFTSUSY 1.9 f8J for the solution of the 
Renormalization Group Equations (RGE) to compute the Supersymmetry mass spectrum at the 
weak scale. 

The two spectrum calculators both consider 2-loop radiative corrections in computing the 
running of masses and couplings between the electroweak and the unification energy scales. They 
differ in the implementation of these corrections, however, and the difference between their results 
reflects the uncertainties on the contributions from higher order radiative corrections. A detailed 
comparison between the RGE calculators can be found in ref. [9], while the uncertainties on the 
resulting predictions of relic density are discussed in ref. [10]. Over most of parameter space, the 
mass spectrum is predicted with an uncertainty of less than 1%, but in the focus-point region 
larger differences are found: in this region of parameter space the results are particularly sensitive 
to the value of the top Yukawa couplings and the uncertainties from higher order corrections are 
significant. 

In Fig. [T]a scan of the (moiiwi/s) plane performed with ISAJET+micrOMEGAs is presented, 
for fixed values of tan/3 = 10, A = 0, and positive fi. A top mass of 175 GeV was used. The dark 
grey region on the left is the scalar tau co-annihilation strip, while that on the right is the focus 
point region with < Qdm- 
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The latter is found at large value of mo > 3 TeV, hence in this scenario the scalar particles 
are very heavy, near or beyond the sensitivity limit of LHC searches. Since mi/2 << mo, the 
gauginos (chargino and neutralino) and gluino states are much lighter. In this scenario the SUSY 
production cross section at the LHC is thus dominated by gaugino and gluino pair production. 

It is instructive to consider the variation of the relic density and the Higgsino mass term /j, along 
a line in the (mo, m 1 / 2 ) plane at fixed m 1 / 2 . This is shown in Fig. [2] for m 1 / 2 — 300 GeV. The left 
plot reports the dependence of fi on mo. When the value of [i drops, the lightest neutralino acquires 
a significant Higgsino component and the relic density decreases (as shown in the right plot) : this 
is the focus point region. The picture shows that ISAJET and SOFTSUSY are in reasonable 
agreement for low values of mo but their predictions diverge as the scalar mass is increased, and 
they find the drop of /i at different values of mo. In addition, the value computed by SOFTSUSY 
never gets low enough to result in an acceptable value of the relic density for this particular choice 
of the parameters. These uncertainties make it difficult to decide whether a given SUSY mass 
spectrum is really consistent with the mass and coupling unification at the high scale assumed 
by mSUGRA. However, this does not prevent us to select a benchmark point with an interesting 
phenomenology, compatible with accelerator and cosmological constraints, and study the ATLAS 
potential to study this model. 

The gluino mass increases with m\/ 2 . It is about 800 GeV for mi/ 2 = 300 GeV, at the bottom 
of the focus point strip allowed by cosmological constraints and accelerator searches. This value 
corresponds to a cross section for gluino pair production at the LHC of about 1 pb. The gluino 
decays to \qq (a chargino or neutralino, and two quarks) followed by the cascade decays of the 
chargino (or neutralino) into the lightest neutralino. These events have the classical mSUGRA 
signature of hard jets and missing energy and can be discriminated from the Standard Model 
background as will be shown in Section [3l 

As one moves upward inside the focus point strip, the gluino mass increases, and it reaches the 
value of 2 TeV for m Y j 2 ~ 900 GeV. This is expected to be roughly the limit of the ATLAS discovery 
potential [11 . The mass of the lightest neutralino and chargino states is much smaller than the 
gluino and squark masses, and the cross section for \X production at the LHC can be larger than 
the cross section for gluino pair production. Since the gauginos are very light, however, their decay 
does not provide a signature with hard jets and missing energy, which would allow discrimination 
against the Standard Model background. The possibility to detect this signal through multi-lepton 
signatures is under study, and is outside the scope of this note. 

No focus point solution is found for tan/3 < 7. For large values of tan/3, the focus point 
solutions move to lower values of mo as shown in Fig. [3] For a fixed value of m\j 2l the gluino 
and gaugino masses are hardly affected. However, the mass of the scalars becomes smaller, and 
the production of qg and qq pairs, followed by the squark cascade decay into gluinos, chargino 
and neutralinos, may be observed at LHC. Note also that the width of the coannihilation strip is 
heavily affected by the value of tan /3. 

The position of the focus point strip is also very sensitive to the value of the top mass. A 
lighter top pushes the focus point strip to smaller values of mo (Fig. |4j. The relation between the 
value of the scalar mass mo for which a focus point solution is found and the top mass is shown in 
Fig. [5] for fixed values of tan j3 = 10, A = 0, fi > and m 1 / 2 = 300 GeV. 

The location of the focus point strip in the (mo,m 1 / 2 ) is much less sensitive to the choice of 
the sign of /i or the value of A, at least in the range between -1000 GeV and 1000 GeV. 

From the considerations above the following point in the parameter space was chosen for the 
detailed study reported in the next sections: 

m = 3550 GeV, m 1/2 = 300 GeV, A = GeV, n > 0, tan/3 = 10 

with the top mass set to 175 GeV and the mass spectrum computed with ISAJET. In Tabled] 
the mass spectrum for this point is given. The scalar partners of Standard Model fermions have 
a mass between 2131 GeV and 3574 GeV. The neutralinos and charginos have masses between 
103.3 GeV and 294.9 GeV. The gluino is the lightest strongly interacting state, with a mass of 
856.6 GeV. The lightest Higgs boson has a mass of 119 GeV, while the other Higgs states have a 
mass at more than 3 TeV, well beyond the expected LHC reach. 

At this point of the mSUGRA parameter space, the total SUSY production cross section at the 
LHC, as computed by HERWIG [T3] at leading order, is 5.00 pb. It is dominated by the production 
of gaugino pairs, x°X° (0.22 pb), x°X ± (3.06 pb), and XX (1.14 pb). 
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Figure 1: The picture shows the regions of the (mo, m 1 / 2 ) mSUGRA plane which have a neutralino 
relic density compatible with cosmological measurements in dark grey. The black regions are 
excluded by LEP. The light grey regions have a neutralino relic density which exceeds cosmological 
measurements. White regions are theoretically excluded. The values of tan/3 = 10, A$ = 0, a 
positive n, and a top mass of 175 GeV were used. The RGE were solved using ISAJET. 



u, m =300 GeV, tan 13=10. m =175, A =0 

^ 1/2 j r— ) top o 





iSAJET +MICROMEGAS 
SOFTSUSYh-MICROMEGAS 



4500 5000 
m„ (GeV) 



1< 2000 2500 3000 3500 4000 4500 5000 

m, (GeV) 



Figure 2: Left plot: dependence of the Higgsino mass term \x on the mSUGRA common scalar 
mass to , for TOiy 2 = 300 GeV, tan/3 = 10, A = 0, a positive fi, and a top mass of 175 GeV. The 
circles are obtained using ISAJET to solve the RGEs, the open squares using SOFTSUSY. Right 
plot: dependence of the neutralino relic density on toq. 
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Figure 3: Same as Fig.[U but for tan/3 = 50 (left plot) and tan/3 = 54 (right plot). 
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Figure 4: Same as Fig.Q] but for a top mass of 172 GeV (left plot) and 178 GeV (right plot). 
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Figure 5: The picture shows the values of the mSUGRA common scalar mass mo and the top mass 
mt which have a neutralino relic density compatible with cosmological measurements. The values 
of mi/2 = 300 GeV, tan/3 = 10, ^4o = 0, and a positive \i were used. The RGE were solved using 
ISAJET. 
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Mass (GeV) 


xS 


103.35 


Sl 


3547.5 


X°2 


160.37 


en 


3547.5 


X°3 


179.76 




3546.3 


xi 


294.90 




3519.6 


xt 


149.42 


h 


3533.7 


xi 


286.81 




3532.3 


9 


856.59 


h 


119.01 




3563.2 


H° 


3529.7 




3574.2 


A 


3506.6 


h 


2924.8 


H± 


3530.6 


h 


3500.6 






ii 


2131.1 






h 


2935.4 







Table 1: The mass spectrum of the benchmark focus point described in the text. 
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Process 


Generator 


Events (10°) 


P r i_t /n — 1 \ 

J i( it(fb ) 


L(cm z s 1) 


SUSY 


HERWIG 


0.15 


30 


10 66 


SUSY 


HERWIG 


1.5 


300 


10 34 


tt 


i\/rr^^ATT o 
IvlCtyJlN-LU 


lb. i 


an 
ZZ 


1 n33 


IF+jets 


ALPGEN 


18.2 


18 


10 33 


Z+jets 


ALPGEN 


6.3 


20 


1Q 33 


bbq 


ALPGEN 


10.9 


0.49 


10 33 


bb + N jets (N > 1) 


ALPGEN 


11.6 


2.4 


1Q 33 


N jets (TV > 2) 


ALPGEN 


2.7 


0.0016 


10 33 



Table 2: Simulated data samples used for the fast simulation studies reported here. 



The production of gluino pairs (0.58 pb) is also significant. The gluino decays into x°qq (29.3%), 
X°g (6.4%), or x ±f lQ (54.3%). The quarks in the final state belong to the third generation in 75.6% 
of the decays. 

The more recent Tevatron data favor a lighter top mass [T^] . As the location of the focus point 
region in parameter space is quite sensitive to the value of the top mass, it is opportune to discuss 
how the phenomenology of our benchmark would change by using a lower value of 172 GeV for 
the top mass. The correct relic density can be obtained lowering the value of mo to 2000 GeV. 
The gluino would then have a slightly lower mass of 818 GeV. The neutralino and chargino masses 
are also slightly smaller, but remain within 10 GeV of the values of those of our benchmark point. 
The leptonic decays of the x° an d X3, which are the basis of the analysis discussed in section [U 
still occur but with a smaller branching ratios (1.9% and 3.3% instead of 3.3% and 3.8%). 

In this scenario, the event rates would be slightly different than in our benchmark scenario, 
but the procedures and the general conclusions of the studies presented in the next sections would 
still hold. 

In addition, the lower value of mo would probably open up the possibility to discover the scalar 
quarks at the LHC . 

3 Inclusive searches 

In this section, the possibility to detect an excess over Standard Model expectations in the produc- 
tion of events with hard jets and large missing energy is investigated. The cuts will be optimized 
to be sensitive to the production of gluino pairs, but they will not rely on any specific decay of 
gluino or its daughters. 

The production of Supersymmetry events at the LHC was simulated using HERWIG 6.55 [13 . 
The top background was produced using MC@NLO 2.31 [2]. The fully inclusive tt production 
was simulated. This is expected to be the dominant Standard Model background for the analysis 
presented here. The Z+jets and the W+jets backgrounds were produced with ALPGEN 2.05 p~5] 
for the hard process, requiring at least 2 jets with px > 40 GeV and \rj\ < 6, and HERWIG for 
subsequent parton shower and hadronization. The vector bosons were forced to decay into leptons. 
The 66+jets background was simulated using ALPGEN for the hard process, requiring \rj\ < 5 and 
Pt > 40 GeV for the b quarks and \r)\ < 3 and pr > 40 GeV for the additional jets (at least one). 
HERWIG was used for the simulation of the subsequent parton shower and hadronization. The 
QCD light jet background is expected to be negligible for the analysis presented here, which require 
either two tagged 6-jets or isolated leptons. In order to verify this assumption, the production of 
events with 3 or more jets with px > 40 GeV was simulated using ALPGEN for the hard process 
and HERWIG for parton shower and hadronization. All jet flavours except 6-jets were included. 

The events were then processed by ATLFAST [16] to simulate the detector response. 

A list of the simulated data used for fast simulation studies is shown in Table [2] The different 
instantaneous luminosities are taken into account through the parametrization used by ATLFAST 
to simulate the detector response, since during operation at low luminosity (10 33 cm _2 s~ 1 ) the de- 
tector performance is not degraded by pile-up effects. The standard parametrization of the ATLAS 
6-tagging performance was used, which assumed a 6-tagging efficiency of 0.6 at low luminosity and 
0.5 at high luminosity, for a rejection factor (the inverse of efficiency) for charm and light quark 
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jets of 10 and 100 respectively. Muons and electrons are considered isolated if they are separated 
by the closest calorimeter cluster in the plane of pseudorapidity and azimuthal angle by at least 
AR = y/ Ai] 2 + A</> 2 > 0.4 and if the transverse energy measured by the calorimeter in a cone 
centered on the lepton and of width Ai? = 0.2 is not larger than 10 GeV (excluding the energy of 
the lepton itself). 

Trigger efficiencies were not taken into account by the detector simulation. The cuts of the 
analyses described in this and in the following sections are more stringent than the selections of the 
ATLAS trigger menu foreseen for operation at the design luminosity of LHC |17j . In particular, 
the cuts on jets and missing energy are more severe than the trigger requirement of at least one 
jet with p t > 70 GeV and Ej Iiss > 70 GeV. Many of the events used in the analysis presented in 
the next section, with two isolated leptons in the final state, would also be selected by the electron 
and muon triggers. A study of the impact of trigger efficiencies on event rates, as well as the 
inclusions of detector effects not described by the parametrized simulation is outside the scope of 
this document. 

The analysis presented in this paper all require a transverse missing energy of at least 70 GeV 
or more, and either two jets with pt > 30 GeV tagged as &-jets or isolated leptons. With these 
cuts, the main contribution from the QCD multi-jet background is expected to come from the 
bb+jets events. This is demonstrated by the results reported in Table [3j which show that after 
these preliminary cuts mentioned above, the light jets production is reduced to less than 10% of 
the heavy flavour background. In the following, the light jet background would not be considered. 



Process 


Events 


^Miss cut 


2 b jets 


1 lepton 


bb+jets 
light jets 


94581.3 
2.69346e+06 


344.296 
193.282 


106.074 
4.8 


20.8 
1.6 



Table 3: Contributions of the multi jet background to the signatures studied in this paper, evaluated 
with ATLFAST events for low luminosity operation. The number of events corresponds to an 
integrated luminosity of 1.6 pb _1 . The third column shows the number of events which have 
^Miss > ^ GeV. The fourth column reports the number of events after the additional requirement 
of 2 tagged b jets, the last column the number of events with Eyr^ , > 70 GeV and one isolated 
electron or muon. 

The most abundant gluino decay modes are g — > x°tt (27.9%), g — > X + tt> (22.0%) and g — > x~ib 
(22.0%B Events with gluino pair production have thus at least four hard jets, and may have many 
more additional jets because of the top hadronic decay modes and the chargino and neutralino 
hadronic decay modes, such as x\ — > XiQQ or xf ~ * XiQQ ■ When both gluinos decay to third 
generation quarks, which happens for 57.2% of the events, at least 4 jets are fe-jets. A missing 
energy signature is provided by the two x.i m f ne final state, and possibly by neutrinos coming 
from the top quark and the gaugino leptonic decay modes. 

These events can be separated from the Standard Model background requiring the presence of 
hard jets and large missing transverse energy. The request of &-jets in the final state suppresses 
the W+jets and Z+jets background so that the dominant surviving Standard Model background 
is the top pair production and the 66+jets production. This request also enhances the signal with 
respect to these backgrounds, as most signal events have four true b-jets (rather than two) in the 
final state. 

The following selections are made: 

• At least one jet with pr > 120 GeV 

• At least four jets with pr > 50 GeV, and at least two of them tagged as &-jets. 

• E T MISS > 100 GeV 

• E llISsl M EFF > 0.12 

1 In the following, whenever the decay g — > tb\~ will be mentioned, the charge conjugate is implied. 
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Sample 


Events 


Basic cuts 


2 6-jets 


M e ff > 1600 GeV 


SUSY (xx) 


a a onn 
44ZUU 


ooy 


on 


U 


SUSY (gg) 


5800 


3522 


1625 


508 


tt 


7.6 • 10 6 


174993 


39816 


506 


TU+jets 


10.1 • 10 6 


62546 


397 


18 


Z+jets 


3.15 • 10 6 


45061 


306 


20 


66+jets 


272 • 10 6 


39579 


12124 


141 



Table 4: Efficiency of the cuts used for the inclusive search, evaluated with ATLFAST events 
for low luminosity operation. The number of events corresponds to an integrated luminosity of 
10 fb _1 . The third column shows the number of events which pass the cuts on jet transverse 
energy, transverse missing energy, and E^ iiss /Meff reported in the text. The number of events 
shown in the fourth column is obtained after the additional requirement that two of the jets are 
tagged as 6-jets. Finally, the number of events which passes all the selection and have an effective 
mass larger than 1600 GeV is shown in the last column. 



Here, the effective mass Meff is defined as the scalar sum of the transverse missing energy and 
the transverse momentum of all the reconstructed hadronic jets. The fraction s — EJ iiss /Meff 
measures the relative importance of the missing energy and the hadronic jet components of Meff- 
In events due to the production of a gluino pair this ratio is on the average larger than in the 
background processes. 

The efficiency of these cuts is reported in Tab. HI The SUSY events are divided in gluino pair 
and gaugino pair production. The selection efficiency for the gg events is 60.7% after all the cuts 
except the request of two 6-tagged jets. This request reduces the selection efficiency to 28.0%. The 
ratio between the number of XX an d gg events is already suppressed by two orders of magnitude 
by the selections on hard jets and missing energy and it becomes negligible after the request of 2 
b— jets. 

The requirement of two 6-jets dramatically reduces the VF+jets and Z+jets backgrounds. Be- 
cause of the large number of true b jets in the SUSY events, it also reduces the signal less than the 
tt and bb backgrounds. These are the dominant sources of Standard Model backgrounds after all 
selections. 

Effective Mass 1 
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• tt, MCatNLO 
□ Z+jetS, ALPGEN 
■ W+jets, ALPGEN 
a bb+jets, ALPGEN 
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Effective Mass (GeV) 



Figure 6: Distribution of the effective mass defined in the text, for SUSY events and the Standard 
Model background, for an integrated luminosity of 10 fb~ . 

The distribution of the effective mass after these selection cuts is reported in Fig. [51 The 
statistics corresponds to an integrated luminosity of 10 fb _1 . The number of signal and background 
events with an effective mass larger than 1600 GeV is comparable and is reported in the last column 
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of Table [H 

From these numbers, it is possible for the focus point scenario under study to estimate the 
minimum integrated luminosity required to observe a deviation from Standard Model expectations. 
The statistical significance of the signal from Supersymmetry is S/y/B = 6.1 for an integrated 
luminosity of 1 fb _1 . These numbers assume nominal detector performances (that is, detector 
commissioning has been completed) and includes only statistical errors. In order to reach this 
sensitivity the systematic error on the background, whose rate is similar to that of the signal after 
all cuts and at large effective mass, will need to be reduced to a level well below a s ^g^(B)/B = 
y/B/B = 12% (for 1 fb" 1 of integrated luminosity). 



4 The di-lepton edge 

As discussed at the end of Section [21 at the selected benchmark point the neutralinos are produced 
either directly or by gluino decays. Despite the lower cross section, the latter mechanism dominates 
after the cuts on missing energy and jets, which are necessary to remove the Standard Model 
backgrounds. The leptonic decays of the second and third neutralino 



xl 



x\i + i~ 



(1) 



x^ 



x\i + i~ 



(2) 



occur with a branching ratio of 3.3% and 3.8% per lepton flavour respectively. The two leptons 
in the final state provide a clear signature. Their invariant mass distribution has a kinematic 
endpoint value equal to the mass difference of the two neutralinos involved in the decay, which is 



m^o = 57.02 GeV m^o 



76.41 GeV 



(3) 



1 / 





Figure 7: Feynman diagrams of the \2 leptonic decay. The decay of the X3 proceeds according to 
a similar diagram. 

The fourth neutralino is heavy enough that the decay x° — > Xi^° is open; the leptons produced 
by this decay have an invariant mass equal to the Z mass and do not allow to measure the neutralino 
masses. 

The leptonic decays of X2 an d X3 proceed according to the two diagrams reported in Fig. [7] 
The resulting distribution of the lepton four momenta is given in Ref. [18] . In the focus point 
region, the diagram with the virtual slepton exchange is negligible, because of the large slepton 
mass. Assuming that only the diagram with the Z exchange contributes, and neglecting the mass 
of the leptons in the final state, we get the following expression [TjJ| for the distribution of the two 
lepton invariant mass: 

£ = Cm ^-^f +^) + (MM) a [_W + m\2M> + (j?) + (,M) 2 ] (4) 
am (m z — m z ) z 

In this formula, C is a normalization constant, /1 = — mi and M = m% + mi, where mi 
and m,2 are the signed mass eigenvalues of the daughter and parent neutralino respectively. At the 
focus point, the mass eigenvalues of the two lightest neutralinos have the same sign, while the \o 
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has the opposite sign. In the decay of the Xo it is thus /i = m{x%)~ m(xi) an d M — m(x2) + m(xi)- 
In the decay of the Xo the role of /i and M is inverted 0. 

The shape of the distribution is different in the two cases, as can be seen from Fig. [8] where the 
distribution of the invariant mass of the two leptons from the two decays is shown. The histogram 
is the true invariant mass of the two leptons and the line is a fit performed with the function of 
eq. H 

The fit provides both the difference and the sum of the neutralino masses. The latter is 
however much less precise, since the dependence of the distribution shape on the lightest neutralino 
mass becomes very weak when this mass is larger than the difference between the masses of the 
neutralinos. The results of the fit to the Monte Carlo distribution are in agreement with the 
true values of the neutralino masses, with the exception of the m^x®) ~ m (Xi) mass difference, 
which is precise enough to be affected by the mass of the muon, which is neglected in eq. HI The 
formula neglects the distortion of the shape caused by the kinematical cuts. The systematic errors 
on the mass difference arising from the analysis cuts on the lepton transverse momentum and 
pseudorapidity have been estimated as follows. When the kinematical cuts of the analysis on the 
lepton momenta (px > 10 GeV and |?7| < 2.5) are imposed, the value of the endpoints obtained 
from the fit of the Monte Carlo Truth distribution changes by 0.20 GeV for the mfyz) — n^iXi) 
difference and by 0.04 GeV for the rn(x.3) ~ rn (Xi) difference. This was taken as the systematic 
error induced on the endpoint by the kinematical cuts. 

The analysis of the simulated data was performed with the following selections: 

• Two isolated leptons with opposite charge and same flavour with px > 10 GeV and \rj\ < 2.5 

• E luss > 80 GeV ' m eff > 1200 GeV, Ej iISS /M EFF > 0.06 

• At least one jet with pt > 80 GeV, at least four jets with px > 60 GeV, and at least six jets 
with p T > 40 GeV 

The efficiency of the various cuts is shown in Table [5] for low-luminosity running conditions and 
an integrated statistics of 30 fb _1 . After all cuts, 411 SUSY and 83 Standard Model events are left 
with a 2-lepton invariant mass smaller than 80 GeV. In 247 of the SUSY events the decay (p} or 
@ is indeed present in the Monte Carlo Truth record. In Fig. [5] the distribution of the invariant 
mass of lepton pairs after all cuts is shown for the signal and the Standard Model background (full 
and dashed lines respectively). 

With these selections, the significance SU SY/ y/SM would be 8.2 with an integrated luminosity 
of lfb -1 (neglecting systematic errors), which makes this channel competitive with the inclusive 
search in terms of discovery reach. 

The background can be estimated from the data using the e + /i~ and /i + e~ pairs. In Fig.[5]the 
distribution of the lepton invariant mass is reported for the same flavour and the opposite flavour 
lepton pairs. Outside the signal region and the Z peak the two histograms are compatible. The 
opposite-flavour distribution was thus used to estimate and subtract the SUSY combinatorial and 
Standard Model backgrounds from the same flavour signal histogram. 

We note that this technique may be used also in the search for a deviation from Standard Model 
prediction. Since part of the SUSY events is removed in the subtraction, the statistical significance 
would be smaller than that obtained with the number of events with opposite sign, same flavour 
leptons. However, the Standard Model contribution is canceled in the subtraction, removing many 
sources of systematics errors on background rates. 

After the flavour subtraction and for a statistics of 30 fb _1 , 179 ±32 SUSY events and 19 ±37 
Standard Model events are left 3 If the invariant mass of the lepton pair is required to be smaller 

2 The second term in parenthesis in Eq. [4] m 2 (2M 2 + fi 2 ), is not invariant for the exchange of fj, and M. This 
term comes from an interference term in the matrix element of the decay and does depend on the relative sign of 
the two neutralino eigenstates. 

3 The quoted error is the uncertainty on the Monte Carlo prediction and it is not equal to the square root of 
the number of events. The expected number of events is computed as b X ((Nsp+ — Ngp_) — (Nqf+ — Nof—)) 
where b is the factor needed to rescale the Monte Carlo statistics to 30 fb _1 , N$f+ and Nsf— are the number of 
same-flavour events with positive weight and negative weight respectively, Ngp+ an d Ngp_ are the corresponding 
numbers for opposite-flavour events. The error on this number is b X y / Ngp+ + Ngp— + Nqf+ + Nof—- For the 
SUSY sample all events have positive weight and b = 1 while for the top sample b = 1.36 and events with both 
positive and negative weights are present. 
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Figure 8: Distribution of the invariant mass of lcpton pairs from the X2 ~ * Xi^ decay (left plot) 
and the X3 ~ * Xi^ decay (right plot). The histogram is the HERWIG Monte Carlo distribution, 
the line is a fit performed with the analytical formula described in the text. 
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Table 5: Efficiency of the cuts used for the reconstruction of the neutralino leptonic decay, evaluated 
with ATLFAST events for low luminosity operation. The number of events corresponds to an 
integrated luminosity of 30 fb~ . The third column contains the number of events after all cuts, 
and the last column reports the number of events with a lepton invariant mass in the signal region. 
SUSY events are divided in gluino pair production with the presence of either the X2 ~~ * Xil + l~ or 
the X3 — * x\l + l~ decay (signal), gluino pair production without these decays (background), and 
the XX production. 
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Figure 9: The full and dashed lines are the distribution of the invariant mass of lepton pairs 
with the same flavour and opposite charge for SUSY events and the Standard Model background 
respectively. The full markers (Supersymmetry) and the empty markers (Standard Model) are the 
distribution of invariant mass of the lepton pairs with opposite flavour and opposite charge. The 
number of events correspond to an integrated luminosity of 30 fb _1 . 
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than 80 GeV, this leaves 170± 25 SUSY events and 8± 13 Standard Model events. The significance 
of the signal from neutralino leptonic decays is 2.5 for 1 fb _1 of integrated luminosity. 

After the flavour subtraction, the Standard Model background is compatible with 0, but still 
contributes to the invariant mass distribution by increasing the statistical fluctuations. The effect 
is small, since even before flavour subtraction the Standard Model contribution is smaller than 
the SUSY combinatorial background. For the high-luminosity studies reported in the rest of this 
section, the Standard Model contribution is not included. 

The flavour subtracted distribution is shown in Fig. [10] for an integrated luminosity of 300 fb _1 
and high luminosity conditions. The presence of two edges is now apparent. The fit was performed 
with the sum of the Xg and x° decay distributions provided by Eq. [H convoluted with a Gaussian 
smearing obtained from the width of the observed Z peak. The fit parameters are the mass of the 
Xi (which is the same for the two decays), the two mass differences x_2 ~ Xi an d X3 ~ Xi> an d the 
normalizations of the two decays. A good fit x 2 is only obtained using the correct values for the 
sign of the neutralino mass eigenstates. 

The values found for the two mass differences are m(x<]) — w(xi) = (57.2 ± 0.4 ± 0.2) GeV and 
m (X3) ~ m (Xi) — (78.1 ± 1.4 ± 0.04) GeV. The first error is the statistical one and the second is 
the systematic error due to the distortion of the distribution arising from the lepton kinematics 
cuts discussed earlier. They are compatible with the true values (eq. [3]). The mass of the lightest 
neutralino is not constrained by the fit, which gives m(xi) = (0.3 ± 2.1) TeV. 

The number of events due to each neutralino decay can be computed as the integral of eq. [5] 
The ratio of the events produced by the two decays can thus be determined from the parame- 
ters measured by the fit. Taking into account the correlation between the parameters, this gives 
N3/N2 — 1.4 ± 0.3. This is compatible from the value N3/N2 = 1.19 which can be computed from 
the gluino, chargino and neutralino branching ratios. 

5 Reconstruction of the gluino decays 

The standard technique used to reconstruct the squark and gluino decays in mSUGRA is the 
combination of the leptons from the neutralino decays with the hardest jets in the event [HI [21]. 
In models for which the two hardest jets in the event correspond to the quarks from the q — > qx° 
decay, it is possible to control the combinatorics from wrong jet associations and reconstruct all 
the masses of the particles in the g — > qq — > X2IQ ~ > U<1<1 — > Xi^QQ decay chain. 

In the focus point region, however, the statistics of lepton pairs is not very high to begin with, 
and the three-body decays of the gluino and the presence of top quarks in most of these decays 
result in a large number of jets (between four and twelve, depending on the gluino and top decay 
modes) which comes from the decay of two gluinos. It is thus very difficult to control the resulting 
jet combinatorics in these events. 

Instead, the strategy used was the explicit reconstruction of a top quark in the event using the 
hadronic decay mode, followed by the selection of tt and tb/tb pairs, with appropriate kinematic 
cuts (in particular, the angular separation) in order to reconstruct the g — > tix° and g — > tbx~ 
decays. The invariant mass distributions of these decays are shown at the parton level in Fig. 1111 

The distribution of the tt invariant mass has four end points, one for each neutralino state. 
Because of the poor experimental resolution on the top reconstructed momentum, they are unlikely 
to be separated from the data. The distribution falls almost linearly to the third endpoint at 
m (5) — 77i(x2) = 696 GeV, and only very few events due to the g — > Xiti decay are found at larger 
invariant mass. 

The distribution of the tb invariant mass has two endpoints at 569.8 GeV and 707.2 GeV, 
corresponding to the difference between the mass of the gluino and that of the two charginos. 
Again, the first endpoint will be very difficult to extract from the data, because of the smearing 
from the finite jet energy resolution. 

5.1 Reconstruction of the g — > x°ti decay 

The reconstruction of the tt invariant mass using the decay tt — > jjbjjb requires the presence 
of six jets, two of which tagged as 6-jets to reduce the combinatorial background. Other jets are 
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Figure 10: Flavour-subtracted distribution of the invariant mass of lepton pairs, for an integrated 
luminosity of 300 fb _ . The fit function is superimposed as a full line; the contribution from the 
X2 and X3 decays are shown separately as a dashed and dotted line respectively. 




Figure 1 1 : Parton-level invariant mass of the tt pairs from the g — > ttx° decays (left) and of the tb 
pairs from the g — ► tb~x~ decays (right). The contributions from the decay to the various neutralino 
and chargino states are shown separately. 
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expected for signal events, from the decay of the other gluino. Thus, events were selected according 
to following cuts: 

• EJ IISS > 120 GeV 

• At least one jet with pr > 150 GeV, at least 8 jets with pt > 40 GeV, at least two of these 
tagged as 6-jets. 

• M E ff> 1200 GeV 

All the pairs of jets which were not 6-tagged, with a transverse momentum pt > 30 GeV, and 
with an invariant mass within ±20 GeV of the nominal W mass were used to build W candidates. 
The combinatorial background is estimated from the events which contain jet pairs in the regions 
A: \m ]0 -{m w -i0GeV)\ < 10 GeV and B: \ mjj - (m w +30GeV)| < 10 GeV. We call them the W 
"sidebands" . The energy and momentum of the jet pairs of each sideband are then rescaled linearly 
by multiplying them by a factor [mw + 2{ m jj ~ (% ± 30GeV))]/rrijj , so that their invariant mass 
lies in the W mass region mw ± 20GeV) [20] . 



j Top invariant mass 
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Figure 12: Distribution of the invariant mass of the top quark candidates (black line) and of the 
combinatorial background estimated from the data, as described in the text (dashed line). The 
plot corresponds to 22 ftV 1 of low- luminosity data and includes both the Supersymmetric signal 
and the Standard Model background. The shaded area indicates the window used to select good 
top candidates. 

Top candidates were found by combining each W candidate with one 6-jet with a transverse 
momentum px > 30 GeV. The invariant mass of top candidates is shown in Fig. [T2l The estimate 
of the combinatorial background is shown as well. The top mass peak is clearly visible over the 
background. The top candidates with an invariant mass within ±20 GeV of the true top mass 
were selected. 

The number of top candidates which passes these selections can be larger than two per event, 
resulting into a large number of possible top pairs combinations. Because of the relativistic boost, 
the top pairs coming from the decay of one gluino have a smaller average angular separation than 
the pairs of top coming from the decay of different gluinos, or those coming from the Standard 
Model ti production. In order to reduce the Standard Model and the SUSY background, the 
angular separation AR (in the (r], <f>) plane) between the two top candidates was required to be 
smaller than 2.5. If more than one top pair passes this selection, we select the combination which 
minimizes 

\/Ki - M t Y + (m t2 - M t ) 2 (5) 

where m tl and m t 2 are the invariant masses of the two top candidates and M t is the true top 
mass. 
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Table 6: Efficiency of the cuts used for the reconstruction of the decay of the gluino into ttx°, 
evaluated with ATLFAST events for low luminosity operation. The number of events corresponds 
to an integrated luminosity of 10 fb~ . The third column contains the number of events which pass 
the inclusive cuts on jets, b-jets, missing energy and effective mass. The fourth column reports the 
number of events with two reconstructed top candidates which satisfy all cuts. SUSY events are 
divided in those with the presence of the g — > decay (signal) , and those without this decay 
(background). 

The number of events which passes the various selections is shown in Table [6] for low-luminosity 
running conditions and an integrated luminosity of 10 fb _1 . The dominant Standard Model back- 
grounds after the inclusive cuts on jets, b-jets, missing energy and effective mass (third column 
of Table O are the tt and the 6&+jets production. The latter is removed when the reconstruction 
of the hadronic decay of two top quarks with AR < 2.5 is required (last column of the table), 
and the dominant background remains the tt production, which is however more than one order 
of magnitude smaller than the signal. 

| Invariant mass best selected tt pair | 




■ SUSY 

n tt 



2 400 600 800 1000 1200 1400 

M INV (GeV) 

Figure 13: Distribution of the invariant mass of the selected pairs of reconstructed top quarks. 
The plot corresponds to 10 fb _1 of low-luminosity data. 

The invariant mass of the top pair is shown in Fig. [T3] for low-luminosity running conditions 
and an integrated luminosity of 10 fb _1 . The statistical significance of the excess of events over the 
Standard Model Contribution is SUSY/VSM = 7.1 for an integrated luminosity of 1 which 
is comparable with the significance expected from the inclusive search and from the di-lepton 
analysis. 

The distribution of tt pair invariant mass for high-luminosity conditions and an integrated 
statistic of 300 fb _1 is shown in Fig. 1141 The high luminosity implies a poorer jet resolution and a b- 
tagging efficiency of 0.5 for the same u-jet mistag probability of 0.01. Only the SUSY contribution is 
included in the analysis for high-luminosity. The contribution from the combinatorial background, 
estimated with top pairs built using the fake W candidates, is also shown. The distribution obtained 
after the subtraction of this contribution is drawn in Fig. 1151 In order to estimate the position of 
the end point, a fit was performed with a polynomial of first order convoluted with a Gaussian. The 
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Figure 14: Distribution of the invariant mass of the selected pairs of reconstructed top quarks. 
The plot correspond to 300 fb _1 of high-luminosity data. 
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Figure 15: Distribution of the invariant mass of the selected pairs of reconstructed top quarks, after 
the subtraction of the background estimated with the side-band technique. The plot correspond 
to 300 fb _1 of high luminosity data. 
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Table 7: Efficiency of the cuts used for the reconstruction of the decay of the gluino into X~tb, 
evaluated with ATLFAST events for low luminosity operation. The number of events corresponds 
to an integrated luminosity of 10 fb _1 . The third column contains the number of events which pass 
the inclusive cuts on jets, b-jets, missing energy and effective mass. The fourth column reports the 
number of events with an accepted tb pair. SUSY events are divided in events with the presence 
of the g — > x tb decay in the Monte Carlo truth (signal), and without it (background). 

Gaussian represents the experimental resolution on the top pair invariant mass, and has been set to 
15% of the endpoint value. The agreement of this function with data is excellent. The endpoint is 
found to be (694 ± 28) GeV, in pretty good agreement with the value of m(g) — m(x2) = 696 GeV. 
The interpretation of the endpoint in terms of this specific mass difference relies on the information 
from the Monte Carlo truth. Still, the endpoint of the distribution of Fig. [T5l provides an estimate 
of the mass difference between the gluino and the neutralino states. 

5.2 Reconstruction of the g — > xtb decay 

The reconstruction of the tb invariant mass using the decay tb — » jjbb requires the presence of 
four jets, two of which tagged as b-jets to reduce the combinatorial background. Other jets are 
expected for signal events, from the decay of the other gluino and/or the chargino. Thus, events 
were selected according to the following cuts: 

• Ej nss > 120 GeV 

• At least one jet with pr > 200 GeV, at least 6 jets with pt > 60 GeV, at least two 6-jets 
with p T > 30 GeV. 

• M E ff> 1200 GeV 

The top candidates were reconstructed as explained in the previous section. They were also 
required to have a transverse momentum px > 150 GeV. If more than one candidate satisfy these 
selections, the one with the closest invariant mass to the top mass was selected. 

Each of the top candidates was combined with the b-jets present in the event. Only the tb 
pairs with an angular separation AR < 2 between the top reconstructed direction and the b-jet 
were considered. In the events where more than one tb pair satisfy these selections, the smallest 
invariant mass of these pairs should be less than the kinematical endpoint of the signal. Since 
combinatorial background pairs with an invariant mass larger than the kinematical endpoint of 
the signal would make the identification of the endpoint difficult, the combination with the smallest 
invariant mass was selected. 

The efficiency of the various cuts is reported in Table[7] For an integrated luminosity of 10 fb — , 
332 SUSY events and 190 Standard Model events pass all the selection cuts, including those on 
the reconstructed tb pair. In 241 of the SUSY events the g — > x tb decay is present in the Monte 
Carlo truth. They are classified as signal in the table, showing the efficiency of the selections for 
these events. 

However, the reconstructed tb pair matches the one from the decay only in 24 of those events. 
In the other SUSY events the reconstructed top does not have a correspondence in the MC truth 
(usually because of the incorrect pairing of jets to build the W candidates), or the tb pair matches 
a true top and bottom quarks originated from a g — > x tt decay or from the decay of two different 
gluinos. The SUSY events thus pass the analysis selections with good efficiency, making this 
channel a promising one for SUSY discovery, but the probability to reconstruct the correct tb pair 
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Figure 16: Distribution of the invariant mass of the selected tb pairs. The plot correspond to 
10 fb _1 of low-luminosity data. 



is relatively small, making it difficult to reconstruct the kinematic endpoint in the invariant mass 
distribution. 

The dominant Standard Model background is again ti. The statistical significance SUSY/VSM 
is 7.6 for an integrated luminosity of 1 fb _1 ; this is comparable to the significance provided by the 
other searches presented here (inclusive, di-lepton and ti analysis). The SUSY/SM ratio is smaller, 
however, than in the leptonic and ti analysis, and the Standard Model contribution to the effective 
mass distribution cannot be neglected. The distribution of effective mass of the tb pairs is reported 
in Fig. [TBI f° r an integrated luminosity of 10 fb _1 . The events are divided in the following classes: 

• SUSY events in which the selected tb pair does indeed correspond to the pair coming from 
the decay g — > x tb. 

• SUSY events, in which the reconstructed top and 6-jet does correspond to top and bot- 
tom quarks in the Monte Carlo truth, but these quarks come from the decay g — > — > 
X°bWbW. These events have not been classified under "SUSY background" because the 
invariant mass of the top with the bottom quark coming from the decay of the other top also 
has four kinematic edges, corresponding to m(g) — m(x ) — m(W). 

• SUSY events in which the 6-jet or the reconstructed top does not have a correspondence in 
the MC truth (usually because of the incorrect pairing of jets to build the W candidates) or 
the corresponding quarks come from the decay chain of two different gluinos. These events 
are classified as "SUSY background". In many of these events the g — > x~tb decay is present 
in the Monte Carlo truth, and are thus classified as signal in table [3 but the reconstructed 
top and bottom do not match those from the decay. The contribution from the fake top 
candidates can be estimated from the sideband distribution reported as a dashed line. 

• Standard Model events 

The contribution from the SUSY and ti backgrounds is large, and the statistics near the ex- 
pected kinematic endpoint is scarce. It is thus not possible to measure the kinematic endpoint. 

Assuming that the contribution of the Standard Model background is reliably estimated, after 
subtraction of the Standard Model and combinatorial backgrounds it may be possible to obtain 
information on the m(g) — rn(x ) difference with larger statistics, but further study is required to 
evaluate this possibility. However, this analysis remains a promising strategy to find evidence of 
an excess of events over the Standard Model contribution. 
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6 Extraction of the MSSM parameters 



In the MSSM, the neutralino and chargino sector of the theory depends on only four parameters: 
the gaugino masses Mi and Af 2 , the Higgsino mass term fi, and the ratio between the Higgs vacuum 
expectation values tan/3. In many models, including mSUGRA, the gaugino masses are unified at 
some high scale. This common value is the parameter rnx/2 m mSUGRA. At the TeV scale the 
relation 

5o' 2 

Mi = ^ T M 2 ~ 0.5M 2 (6) 

holds, where g and g are the electroweak coupling constants. This relation reduces the number 
of free parameters to three unknowns. 

The constraints placed by the dilepton edge analysis on the neutralino mass spectrum can be 
used to determine the values of Mi, \x and tan/3 which are compatible with the experimental 
measurements. 

The ISAJET 7.71 code [7] was used to compute the masses and branching ratios of Supersym- 
metric particles as a function of the soft Supersymmetry breaking parameters. This provides, in 
particular, the values of neutralino and chargino masses as a function of Mi, M 2 , fi and tan/3. 

For each point of parameter space the mass differences AM 2 = m(x 2 ) — m (Xi) an d AM3 = 
TO (X;j) — m (Xi) were compared to the measured values AM 2 ex ^ = 57.2 GeV and AM 3 exp = 
78.1 GeV obtained using 300 fb _1 of simulated data, as discussed in Section [4] The following % 2 
was evaluated: 



X 2 = (AM 2 - AA/ 2 xp ) 2 /ct 2 + (AM 3 - AM 3 exp ) 2 /cj 2 - 2r(AM 2 - AM 2 P )(AM 3 - AM 3 exp )/cr 2 CT 3 

(7) 

where er 2 = 0.4 GeV and 03 = 1.4 GeV are the errors on the measurements of AM 2 and AM3 
respectively, and r = 0.038 is the correlation coefficient. 

The points in the parameter space which have a x 2 probability larger than 0.05 were selected. 
In addition, the same sign was required for the x 2 and Xi mass eigenstates and the opposite sign 
for the X3 an d Xi mass eigenstates, since a good fit of the di-lepton invariant mass distribution 
can be obtained only under these hypothesis, as discussed in Section |U 

A scan of the values of the parameters Mi, fj, and tan/3 was then performed, using relation[6]to 
fix the value of M 2 , to find the parameter space which gives masses compatible with the simulated 
experimental constraints determined in section [5] The scan was performed in the range 

45GeV < Mi < 500GeV, 45GeV < |u| < lOOOGeV, 1 < tan/3 < 65 (8) 

in steps of 2 GeV for Mi and /i, for 16 values of tan/3, and for both positive and negative values 
of fx, so that a total of 3.4 x 10 6 values in the parameter space were considered. The region in this 
space which is compatible with the simulated ATLAS constraints is shown in Fig. [T7] for positive 
values of u and in Fig. [18] for negative values of \i. 

The upper left plot shows the allowed values of Mi and /i. Each parameter can vary in a 
relatively large interval: 97 GeV < Mi < 136 GeV and 138 GeV < u < 182 GeV for positive 
values of fj,, 116 GeV < M x < 157 GeV and -208 GeV < \i < -168 GeV for negative values of fj,. 
However, the ratio is more constrained to the range 1.32 < fi/Mi < 1.46. The value of tan/3 is not 
constrained by the experimental data. 

The value of the ratio between [i and Mi is interesting since in most of mSUGRA space 
u >> Af 2 ~ 1M\, With this hierarchy between the parameters, the lightest neutralino is almost a 
pure Bino. The focus point region of mSUGRA, instead, is characterized by u ~ Mi [5]. In this 
case the mixing between gauge eigenstates is nearly maximal and each mass eigenstate receives a 
significant contribution from all the gauge eigenstates. The Higgsino component of the Xi allows 
rapid s-channel neutralino annihilation, which is the mechanism which reduces the relic density in 
the early universe. The measurement of the two leptonic edges allows to establish this scenario. 

For any given set of values of the parameters Mi , M 2 , /i, and tan /3 the masses of the neutralinos 
and their gaugino mixing angles can be computed, so these quantities are also constrained by the 
data. For the values of the parameters which are compatible with the dilepton invariant mass data 
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Figure 17: Values of the MSSM parameters fx, Mi and tan/3 which are compatible with the 
experimental constraints on the neutralino mass spectrum obtained with 300 fb -1 of data and at 
95% confidence level, for positive values of \i. 
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Figure 18: Values of the MSSM parameters fx, Mi and tan/3 which are compatible with the 
experimental constraints on the neutralino mass spectrum obtained with 300 fb _1 of data and at 
95% confidence level, for negative values of /i. 
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at 95% C.L. the mass of the lightest neutralino lies in the range 67 GeV < m(xi) < 156 GeV. 
For any given value of the mass of xli the mass of the next two lightest neutralinos is precisely 
constrained from the measurements of the two edges of the dilepton distribution. 

7 Conclusions 

A study of the ATLAS potential of detecting and measuring Supersymmetry in the focus point 
scenario has been presented. For the selected point in the parameter space the observation of 
an excess of events over the Standard Model expectations should be observed rather easily; the 
time needed for discovery would probably be determined by the understanding of the systematics 
related to the detector response and the knowledge of the Standard Model backgrounds rather than 
the statistical significance. Several channels can contribute to the discovery of an excess of SUSY 
events with a comparable statistical significance: with appropriate kinematic cuts, the presence 
of Supersymmetry physics may become manifest through an excess of events with hard jets, large 
missing energy and 6-jets (section [3]), events with hard jets, large missing energy and opposite-sign 
electron or muon pairs (section [4j, and events with hard jets, large missing energy, a top and a 
bottom quark or two top quarks (section[5]). In each of these channels, the contribution from SUSY 
events has a statistical significance between 6.1 and 8.2 standard deviations with 1 fb _1 of data. 

In the focus point region, the neutralino leptonic decays proceed through a direct three-body 
decay in which the virtual slepton exchange is negligible, because of the large scalar mass. An 
analytical expression was derived in section U for the distribution of the invariant mass of the 
resulting lepton pairs. This formula was used to fit the distribution of opposite-sign lepton pairs 
obtained with 300 fb^ 1 of simulated LHC data. Two kinematic edges, measuring the m(xV)~ m iXi) 
and m(x2) ~ m ix.i) mass differences, are measured with a precision of the order of 1 GeV. 

The constraints that this measurement would place on the MSSM gaugino sector parameters 
have been discussed. These constraints are such that a scenario with a large mixing in the neutralino 
sector and a relic neutralino density of the same order of magnitude as the Dark Matter abundance 
would emerge from the LHC measurements. 

The gluino decay into ttx° and tbx + can be studied through the reconstruction of the ft and 
tb invariant mass. While a precise reconstruction of the corresponding kinematic edges does not 
seem to be possible, evidence for these decays and an estimate of the mass difference between the 
gluino and gaugino states may be extracted from the experimental data. 



21 



Acknowledgments 



The authors wish to thank G. Polesello (INFN Pavia) for his helpful suggestions on the analysis 
and on the preparation of this document. They also thank D. Tovey (Sheffield), F. Paige (BNL) 
and I. Hinchliffe (LBNL), as well as M.M. Nojiri (KEK, Tsukuba), for useful discussions and 
advices. 

The authors are also grateful to S. Asai and the colleagues of the University of Tokyo for pro- 
viding the ALPGEN 4-vector files used in the production of the W+jets and Z+jets backgrounds, 
and to B. Gjelsten (Bern) for the production of the Z+jets fast simulation ntuples. 

This work has been performed within the ATLAS collaboration making use of the detector fast 
simulation and tools which are the result of collaboration-wide efforts. 

References 

[1] See H.P. Miles, Phys. Rev. 110 (1984) 1 and references therein. 

[2] See the LEP (ALEPH+DE LPHI+L3+OPAL experiments) SUSY working group page, 
http://lepsusy.web.cern.ch/lepsusy[ and references therein. 



[3] D.N. Spergel et al., astro-ph/0603449 



C. L. Bennett et al., Astrophys. J. Suppl. 148 (2003) 1. 

D. N.Spergel et al., Astrophys. J. Suppl. 148 (2003) 175. 

[4] J.R. Ellis et al., Phys. Lett. B565 (2003) 176. 

H.W. Baer and C. Balazs, J. Cosmol. Astr, Phys. 0305 (2003) 006. 



[5] J.L. Feng, K.T. Matchev, and T. Moroi, Phys. Rev. Lett. 84 (2000) 2322 |hep-ph/9908309 . 
J.L. Feng, K.T. Matchev, and T. Moroi, Phys. Rev. D61 (2000) 07500 5 [hep-ph/9909334 
J.L. Feng, K.T. Matchev, and F. Wilczek, Phys.Lett.B482 (2000) 388 |hep-ph/0004043|. 
J.L. Feng, K.T. Matchev, and F. Wilczek, Phys. Rev. D63 (2001) 045024 |hep-ph/0008115] . 

[6] G. B. Bclanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 174 
(2006) 577 |hep-ph/0405253| . 



[7] H. Baer, F.E. Paige, S.D. Prototopescu and X. Tata, hep-ph/0312045 



[8] B.C. Allanach, Comput. Phys. Commun. 143 (2002) 305-331, hep-ph/0104145 



[9] B.C. Allanach et al., JHEP 0303 (2003) 016 hep-ph/0302102 and references therein 



[10] B.C. Allanach et al.,|hep-ph/0402161| 



G. Belanger, S. Kraml, and A. Pukhov, Phys. Rev. D72 (2005) 015003 hep-ph/0502079|. 



[11] ATLAS Physics Performance and Detector Technical Design Report, CERN/LHCC 99-14 
(1999). 



[12] Tevatron electroweak working group, hep-ex/0608032 



22 



[13] G. Corcella et al., JHEP 0101 (2001) 010. 
S. Moretti et al., JHEP 0204 (2002) 028. 
G. Corcella et al., |hep-ph/0201021p . 



[14] S. Frixione and B.R. Webber, JHEP 0206 (2002) 029 [hep-ph/0204244]. 

S. Frixione, P. Nason and B.R. Webber, JHEP 0308 (2003) 007 [hep-ph/0305252| . 



[15] MX. Mangano, M. Moretti, F. Piccinini, R. Pittau, A. Polosa, JHEP 0307 (2003) 001 
[hep~ph/0206293| . 

M.L. Mangano, M. Moretti, R. Pittau, Nucl.Ph.ys. B632 (2002) 343 |hep-ph/0108069] . 
F. Caravaglios, M. L. Mangano, M. Moretti, R. Pittau, Nucl.Phys.B539 (1999) 215 
hep-ph/9807570] . 

[16] E. Richter-Was, D. Froidevaux and L. Poggioli, ATL-PHYS-96-079 and ATL-PHYS-98-131. 



[17] ATLAS High-Level Trigger, Data Acquisition and Controls Technical Design Report, 
CERN/LHCC 2003-022 (2003). 

[18] M.M. Nojiri and Y. Yamada, Phys. Rev. D60 (1999) 015006. 

[19] S. Montesano, "Ricerca di particelle supersimmetriche nell'ambito dell'esperimento ATLAS", 
diploma thesis, University of Milano (2006). 

[20] J. Hisano, K. Kawagoe and M. Nojiri, Phys. Rev. D68 (2003) 035007. 



[21] C.K. Gjelsten et al., ATL-PHYS-2004-007. 
C.K. Gjelsten et al., JHEP 12 (2004) 003. 



23 



